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ЛЕ ЧАН ТХАНГ 


СИНТЕЗ НЕЛИНЕЙНЫХ РЕГУЛЯТОРОВ ПЕРЕМЕННОЙ 
СТРУКТУРЫ ДЛЯ ОДНОГО КЛАССА НЕЛИНЕИНЫХ ОБЪЕКТОВ 


В данной работе основное внимание уделяется проблеме построения управления в 
критической области, где нельзя определять управление так же, как в остальных 
некритических областях. В целом в работе получен интересный теоретический ре- 
зультат, который позволяет предложить новый тип регуляторов для нелинейных 
объектов управления. Этот регулятор называется нелинейным регулятором пере- 
менной структуры. 
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Введение. На практике часто встречаются нелинейные объекты управле- 
ния. Синтез нелинейных регуляторов для таких объектов всегда затрудни- 
телен, так как отсутствует достаточно эффективный метод в общем случае 
[1, 2]. Для решения этой проблемы нами предложен метод синтеза нелиней- 
ных регуляторов с переменной структурой закона управления [3, 4]. 

В теории управления известны регуляторы переменной структуры, 
в которых коммутируются линейные блоки [5]. Эти регуляторы строятся на 
основе фазовых портретов блоков, поэтому в них обычно коммутируются 
линейные блоки второго порядка. Это обусловливает узкие возможности 
подобных регуляторов, в особенности для синтеза систем управления нели- 
нейными объектами. 

Предложенный в [3, 4] метод синтеза характеризуется тем, что 
коммутируются нелинейные законы управления. Поэтому соответствующие 
регуляторы называются нелинейными регуляторами переменной структуры 
(НРПС). Синтез коммутируемых нелинейных законов управления ведется на 
основе функции Ляпунова, что значительно расширяет возможности управ- 
ления нелинейными объектами. Алгоритм НРПС является достаточно слож- 
ным, но он легко реализуется на промышленных микроконтроллерах. Та- 
ким образом, метод, предложенный в [3, 4], позволяет использовать широ- 
кие возможность современной вычислительной техники. 

Постановка задачи: построение алгоритма нелинейных регулято- 
ров переменной структуры в критической области. Рассмотрим 
объект управления, который описывается нелинейной системой дифференци- 
альных уравнений 

я= ЕОьи,), (1) 
где хе А” - вектор состояния; №, - скалярное управляющее воздей- 


ствие; А’(х,и,) - нелинейная дифференцируемая вектор-функ- 
ция, (0, 0)=0. 

Задача синтеза нелинейного регулятора для объекта (1) заключает- 
ся в определении управления # ; = и, (х%), т.е. в определении такой нели- 
нейной обратной связи, чтобы положение равновесия х= О являлось 
устойчивым в смысле Ляпунова в целом [2]. 

Решение. Предполагается, что нелинейность системы (1) удовлетворя- 
ет условию 
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9'Е(х,и ,) 
и =О при 122, (2) 
9 иг 
тогда эту систему можно представить [3] следующим образом: 
х= Ах+ Л) Боди, , (3) 


где „4, - постоянная матрица линейного приближения вектор-функции 
Е’(х, 9) в точке х= 0; Л(х) - некоторая нелинейная вектор- 
функция. Подчеркнем, что уравнение (3) является точным представле- 
нием уравнения (1) при условии (2). 

Если „4 имеет собственное число с положительной вещественной 
частью, то полагаем И, = И + и,, где м, 1, - новые управления. Предпо- 
ложим, что пара 4, ,ЬБ(х) вполне управляема при всех хе В” и что су- 
ществует управление и, = К(х)х такое, что постоянная матрица 
А= А + 6х) К (х) имеет все собственные числа с отрицательными веще- 
ственными частями. При таком управлении 2, система (3) принимает вид 

х= Ах+ Г(х)+ Б(хи. (4) 

Изложенное позволяет считать, что постоянная матрица „4 обес- 
печивает асимптотическую устойчивость линейной части (4). В этом случае 
существует положительно определенная симметрическая матрица Р, ко- 
торая является решением уравнения Ляпунова: А7’Р+ РА=-С, где С 
— положительно определенная матрица [2]. Возьмем в качестве кандидата 
в функции Ляпунова для системы (4) положительно определенную квадра- 
тичную форму 

И) =х'Рх (5) 
и найдем её производную по времени вдоль траекторий системы (4): 
Г) =х'Рх+ х'Рх= 
= -х'Сх+ 2[х' РУ(х) + х'РЬ(хуы]. (6) 

В соответствии с известной теоремой Ляпунова [2, 5], для обеспе- 
чения асимптотической устойчивости положения равновесия х = 0 систе- 
мы (4) необходимо выбрать управление м = м(х) так, чтобы И(х) < 0 
на траекториях этой системы. 

Пусть &, &, - малые положительные числа. Тогда, в соответствии 


с работами [3, 4], алгоритм НРПС включает следующие частичные законы 
управления: 


1. Если при некоторых хе К”, ||х ||< &, имеет место неравенство 
х'Р/(х)< ах'Сх и ЧЕ [0, 0,5), то и(х) = О, т.е. 

и(х) = О при х’Р/(х) < ах'Сх или |х|к&,. (7) 
При этом согласно (6) 

И(х) = -х'Сх+ 2х! Р/(х) < -х'Сх+ 24а х'Сх<0, 

т.е. является отрицательно определенной функцией. 

2. Если х’РХ(х)> ах’Сх, (Е [0, 0,5), || х|> Е, и выполня- 
ется условие | х’РЬ(х) |> =, то 
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и(х) = - х'Р/(х)/х'РЬ(х). (8) 
В этом случае согласно (6) производная Го) =-х'С, х <0 так же. 


Если же | х’РЬ(х) |< = , то управление, определённое по (8), 


может принимать недопустимо большие значения. Поэтому необходимо 
найти новое управление в этой области. 
С учетом соотношения (7) предлагается следующее определение 


критической области: область „Е №” называется критической, если 
при Х0 кр Выполняются условия: 
Их =,, х^Р/(о) > в х'Сх, |х'РЫЮ ке. (9) 
Перейдем к формированию управления в критической области. 


В дальнейшем предполагается цифровая реализация НРПС, поэтому 
далее все непрерывные величины заменяются их дискретными значения- 


ми: 8 Ц, х(Ю > х,, и(х) > и,. При этом период квантования по вре- 
мени АЕ предполагается достаточно малым, так что 
х(В) = х(КА В = ХЕ) = х, при всех 1Е [^АЕ, (+ ПАЙ; 
ЖЕ 5 
Пусть в момент времени 2,„_, траектория системы (4) пришла в 
точку х„_: , которая находится вне области © к» (9). В этой точке, со- 
гласно (8), управление 2,„_, определяется по формуле 
и (10) 
Предположим также, что в следующий момент времени 
11 =, +А ЕЁ траектория системы (4) под действием управления (10) попа- 


дает внутрь области © кр. 

Будем искать управление И, такое, чтобы производная 
ГС, ) < О, т.е. такое, при котором значение функции Ляпунова в следую- 
щий момент времени 2„,, было меньше, чем её значение в момент време- 
ни Е, . Другими словами, при искомом управлении И, должно выполнять- 
ся неравенство 


Рена РА (11) 
Учитывая малое значение /\2 , примем, что 
Хы = х, + Х,А &, где х, = Ах, + Ё(х,)+ Б(х, и. (12) 


Поставим >Х,„., из первого равенства (12) в (11). В результате будем 
иметь 
)РАЛОРО ЕАО Рх, 
ь хРх+х, РхА 1+ Ах. Рх. + Ах. Рх.А1< х,.Рх, 
& 2х Рх+АЁХ,Рх, <0. 
Сократив обе части последнего неравенства на /\ 2 , получим 
2х. Рх,+Аах,Рх, <0. (13) 
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Далее подставим х,„ из второго равенства (12) в (13). В результате полу- 
чим (для большей наглядности выкладок Б(х, ) далее заменено на Б) 
у 
2х, Р[Ах,+ Г(х,)+ Би, ]+ 

+А Ах, + Г(х,)+ Би, Т Р[Ах, + Г(х,)+ Би, ]< 0. 
Это неравенство эквивалентно следующему неравенству второго порядка по 
и, : 

АЕ" РЬ(и,)? + 24А ЦАх, + Г(х, Г РЬ+ х’РЬи, + 

+ 2х + АЩАХ, + Г(х, А РГАх, + Рок 0. — (14) 
Для того чтобы существовало управление И, , необходимо и достаточ- 


но, чтобы неравенство (14) имело вещественное решение. Другими словами, 
его дискриминант А (х) должен быть положительным, т.е. необходимо, чтобы 


А (х)> 0, 
Учитывая структуру полинома в (14), имеем: 
А (,) = {{[Ах,+ ЛОТ РЬУ - 
- ГРЫ Аж, + Г, РГАх, + Го) 2+ 
+24 [Ах, + ОГ РЬ х. РВ - 
-[Ах, + До УГ РГ РЬА 1+ [х/РЬ? > 0. 
Отсюда следует, что искомое управление всегда существует при достаточ- 


но малом /\2. 
При выполнении условия (15) из уравнения (14) находим 


- вы + Ло] РЬ+ жРЬ А) 


АСЬ'РЬ 


Искомое управление 2, находится на интервале (и, и.) ‚ Т.е. 


(15) 


аа 
Можно, например, полагать, 
- А ЦАх, + Го, РЬ+ хГРЫ- 1х, ) 
—- АБ РЬ о НСО 


Таким образом, алгоритм выбора нелинейного управления (>) 
объектом (4) имеет следующий вид: 








Опри  х’Ю/)х ‹ахСхили || ||х <, 
- х'Р/(х)ИвИРЬ ()хРГх „ИК Схих| &| хРЬ Е 
и(х)= - В й = 17 
{А НАх+ ГО) — хтРЫ- у /А © ре (17) 
АЕБ’ РЬ 
х' РА ‹ах'Сх, хи = | |хРЬ 
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Найденное управление (17) решает поставленную задачу, так как 
положение равновесия х= 0 системы (4), (17) является устойчивым по 
Ляпунову в целом, в силу теоремы 4.1 из работы [2. С. 20]. 

Пример. В качестве примера синтеза НРПС возьмем задачу устра- 
нения хаотических колебаний в управляемой системе Рёсслера [1]. Указан- 
ная система Рёсслера определяется уравнениями: 

м =-Х, - %,, 
ЖЕ, (18) 
ЕАО ни 


В данном случае @ = 0,15; В = 0,2; с= 10. Если представить систему (18) 
в виде (3), то соответствующие матрица, нелинейность и вектор будут иметь 
вид: 


О Саи 1 0 0 
А= 1 0,15 0 л ;1(х)= 0 ;5-= 0; 
00 пл -10 а 1 
при и, = но 2щи, (19) 


Матрица „4, не обеспечивает асимптотическую устойчивость линейной 
части (18), поскольку имеет собственное число с положительной веществен- 
ной частью. Пара (.4, , 6) вполне управляема. Возьмем управление 


и (х)= 20,9725х + 12,4959х, + 1,85х. , (20) 
при котором матрица = 4, + Би, (х) из системы (4) принимает вид 
0 -1 -1 
= 1 0,15 0 


20,9725 12,4959 -8,15 


Как будет показано ниже, положение равновесия нелинейной си- 
стемы (19), (20) является неустойчивым, поэтому найдем стабилизирующее 
нелинейное управление м = и(х) . Переходя к его синтезу, примем значе- 


ния параметров & = 0,3, &=0)1, &, =10°, 1= 0,99, ДЕ= 0,001 и 
обозначим: 
НС®) = х'Р Г(%) = (-0,1359х, + 0,0034х», + 01394х. )х. хз, 
ЕО СОЗ а), 
[(х) = х'РЬ= - 0,1359х + 0,0034х, + 1,1394х., 
А (х) -— определяется по формуле (15). 
Тогда по формуле (17), будем иметь 
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0 при Н(5)‹ Цх) или |2; 
- Вх) при Вх)? Цх, Ри [10)} Е; 


и(х): | - | Ах ГО РВ 1 099.1 () 
| ЕБТРЬ —- 
при Н(х)? Ц»), [хи | КЕ. 


Моделирование системы, заданной уравнениями (19), (20), (21), 
проводилось в среде МаНаБ. На рис.1-3 показаны результаты моделирования 


при начальных условиях х, = [- 16,- 50, 70]. На рис.1 представлены гра- 








фики переменных системы (19) только с управлением 14, (20). Как видно, 
переменные х1, хз неограниченно возрастают, т.е. положение равнове- 
сия х = 0 системы, заданной уравнениями (19), (20), является неустойчи- 
ВЫМ. 

На рис. 2 представлены графики переменных системы (17) с управле- 
ниями 2, (20) и м(х) (21). График соответствующего управления м(.х) 
показан на рис. З,а. В этом случае, согласно рис. 2, переменные системы за- 
Тухают, т.е. положение равновесия х = О система, заданная уравнениями 
(19), (20) и (21), является устойчивым по Ляпунову в целом. 


х10“ х10”° 


15 20 





10 0 

















0 Е 0.5 


Рис. 1. Изменение переменных состояния системы, 
заданной уравнениями (19), (20) 
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Рис. 2. Изменение переменных состояния системы, 
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заданной уравнениями (19), (20) и (21) 
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Рис. 3. Управление (21) системы, 
заданной уравнениями (19), (20) 


На рис. 3,6 показан график управления в увеличенном масштабе в 
интервале времени 0,5-0,7 с. Подчеркнем, что согласно рис. 2, несмотря 
на колебания управления, переменные состояния изменяются плавно. 

На рис. 4, а дан портрет системы (19), (20), (21) в трехмерном про- 
странстве, а на рис. 4,6 - 4,г- проекции фазовых траекторий на плоскости 
ХХ; ХХ, Хз и Хх, , х. соответственно. 


10 т т т 
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Рис. 4. Фазовый портрет и проекции 


Из графиков видно, что управление (21) является ограниченным, а 

положение равновесия системы, заданной уравнениями (19), (20), (21), - 
устойчивым по Ляпунову. 
Выводы. Предложенный алгоритм работы нелинейного регулятора пере- 
менной структуры обеспечивает эффективное управление нелинейными 
объектами, и может быть распространён на случай объектов с несколькими 
управлениями. 
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